THE LOWEST-DEGREE POLYNOMIALS WITH 
NON-NEGATIVE COEFFICIENTS 



TOMAS KEPKA AND MIROSLAV KORBELAR 

Abstract. A polynomial p 6 R[as] is a divisor of some polynomial 
0^/6 R[i] with non-negative coefficients if and only if p does not have 
a positive real root. The lowest possible degree of such / for a given p is 
known for quadratic polynomials. We provide it for cubic polynomials 
and improve known bounds of this value for a general polynomial. 



1. Introduction 

Polynomials with non-negative coefficients appear naturally in the ana- 
lysis of stochastic context-free grammars (with applications to natural lan- 
guage processing [3 [9]), control theory (testing of stable polynomials), opti- 
mization and semialgebraic geometry (application of Polya's theorem) and 
algorithmic game theory (bounding the price of anarchy). They are also of 
interest in approximation theory [131 G2] , rewriting systems [17^ YM\ , num- 
ber theory [21 US] or graph theory [20] and have application in biochemistry 
[H [21] and electronics [16] . We will mention here yet another natural mo- 
tivation for studying of these polynomials - the ordered domains and their 
connection with semirings. 

A commutative and unitary domain R with a non-identical (partial) order 
■< on R is called ordered if (1) a<b^>a + c<b + c and (2) (a < b & ■< 
c) =^ ac H be for all a,b,c £ R. In real algebraic geometry the ordered 
domains with additional condition ^ a 2 for all a £ R are studied. This 
property is a weaker form of the linear ordering, but still there are a lot of 
rings (e.g. algebraically closed fields), where the last mentioned condition 
can not be fulfilled, although the definition using only (1) and (2) provides 
many non-linear orders generally. 

In an ordered ring R (necessary of characteristic 0) is the order ^ uniquely 
determined by the set P(^) = {a £ R | -< a} of positive elements. 
This set is a subsemiring of R (i.e. non-empty, closed under addition and 
multiplication) and does not contain zero. For a £ R there exists an order 

on R such that -< a if and only if the semiring P = {/(a)|0 ^ f £ 
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x ■ No[x]} does not contain 0. Since char(i?) = 0, we can consider the 
minimal polynomial ^ minQ(a) G Z[x] of a, if it exists (otherwise put 
minQ(a) = 0). Then we get the following characterization: 

1.1. Theorem. Let R be a commutative domain and a G R. Then there 
exists at least one order X on R such that -< a if and only if the polynomial 
minq)(a) has a real positive root. 

This theorem is an immediate consequence of the following assertion: 

1.2. Theorem. [141 fTI~l [15] Let p G R[x] be a polynomial with no positive 
real roots. Then there is ^ h G Q[x] such that the polynomial f = hp has 
all its coefficients non-negative. 

Probably the first known proof of 11.21 appeared implicitly in papers of 
Poincare |14j and Meissner [IT]. Moreover, by Polya's theorem [15] one can 
choose h = ±(x + l) k with k G N sufficiently large in ll.21 

Apart of this, for that a G R which does not allow the appropriate or- 
dering, the structure of the Q + [x]-semimodule M = {/ G Q + [x]\f(a) = 
0} Q Q. + [x] can be studied (here F + is the semiring of usual non- negative 
numbers of a field F C R and F + [x] is the semiring of polynomials over x 
with coefficients from F + ). In fact, this is a particular version of a general 
problem: 

Having a field F C R and an ideal L C ¥[x], determine the ¥ + [x]- 
semimodule L + := L D ¥ + [x]. 

In this paper we recall and combine known results in connection with this 
problem, especially concerning the estimations of the minimal degree of all 
non-zero elements in L + , and make an improvement of them. It is worth to 
mention that despite the wide range of applications, the polynomials with 
non-negative coefficients are not studied systematically and many papers 
concerning them were done by authors working in different fields of math- 
ematics and using various terminologies. This probably is the reason why 
the results do not follow up each other, miss citations and several of them 
were "rediscovered" repeatedly (e.g. Theorem ll.2p . 

2. Polynomials with non-negative coefficients 

Through this paper, let for a semiring S and an S-semimodule M be (X) s 
the subsemimodule of M generated by a set X C M. 

Let F be a subfield of R, p G ¥[x] and / = pF[x]. First recall, that 
for a given n G No the set M n = {/ G I + |deg(/) < n + deg(p)} is an 
intersection of the vector space (p, xp, . . . , x n p) ¥ over F with the convex F + - 
cone {/ G F + [x]|deg(/) < n + deg(p)}. Using basic knowledge of convex 
geometry one without difficulty gets that M n is thus generated by a finite set 
of " extremal rays" or more precisely, there is the smallest (up to F + -multiples 
of elements) finite generating set X n C F + [x] of the F + -semimodule M n . 
These generators can be computed via linear equations using only rational 
numbers and coefficients of p and therefore they do not depend on the choice 
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of the extension field. Similarly, there is the least generating set X (up to 
F + -multiples of elements) of the F + [:r]-semimodule / + {X may be obtained 
from U n6 N ^n by letting out those polynomials that are x fc -multiples (for 
k G N) of other polynomials in U n ^ X n ). It is also not difficult to show that 
I + is finitely generated as an F + [x]-semimodule if and only / is a monomial 
ideal. 

For a chosen polynomial it is therefore possible (more or less efficiently) to 
compute generators of / + up to any given degree. But to determine explicitly 
I + in general seems to be a difficult task, since even the estimations of the 
least possible degree of a non-zero polynomial in I + are known only very 
roughly. 

2.1. Definition. For O^pG denote deg + (p) the least possible degree 
of a non-zero polynomial in R + [x] divisible by p (if such exists), otherwise 
put deg + (p) = oo. Further put D(p) = deg + (p) — deg(p) and V + (M.) = 
{0}U{/GM[x]|(VaG (0,oo)) /(a) > 0}. 

Polynomials of minimal degree in (p]R[2;]) + for p being quadratic were 
determined by Danes PJ. His proof used Minkowski- Farkas theorem for 
systems of linear inequalities. This result was later proved by Motzkin and 
Strauss [12] by a simpler argument using cones of complex numbers. Finally, 
Harnos [7] described completely the generating set of the R + [x]-semimodule 
(pR[x]) + for p G R[x] \ M + [x] being quadratic. 

2.2. Theorem. [21 [7] Let p = x 2 + ax + b = (x — ui)(x — ul) G M[x], where 
U) = r(cos if + i sin <p), r > and if € (0, it). Then 



deg,(p) 



"7T 









7T 



arcsin( A /l - |^) 



Moreover, there is f G (pIR[x]) + with all coefficients positive such that 
deg(f) = deg + (p) + 1. 

We derive the value of deg + for cubic polynomials. Let us recall that 
every polynomial p of degree d > 1 has a well known correspondence to 
a linear homogeneous recurrence relation with constant coefficients. All 
sequences that are solution of this relation form a vector space of dimension 
d with a particular basis (a^)^ , . . . , (a\ )£^o ma de up with help of roots 
of p and derivations. Now put (vfc(p) =)vjt = (a^ 1 , . . . , a£ ) G M. d for 
k G No- The vectors Vo,. . . , v^_i are linearly independent and thus f(-) = 
det(vo, . . . , v^_2, • ) is a non-zero functional. Put nT(p) := minjfc G {d, d + 
1, . . . }|y(vfc_i)</?(vjfc) < 0} G N U {oo} and n^(p) := rv(ft), where p{x) = 
x p(h)- Since there is ^ f G (pM[x]) + of degree m if and only if v m G 
(vo, . . . , v m _i) R+ and since (p corresponds to a particular halfspace, it follows 
that n^(p) < deg + (p). Further, as one easily checks, the number n^(p) will 
be the same if we choose any other basis than (a- )f2. , ■ ■ ■ , (aj (e.g. 
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such one where the initial values correspond to the identity matrix). We 
summarize now our conclusions into 12.31 and I2.4L 

2.3. Proposition. Letp G be of degree > 2. Then max{n^(p), n^(p)} < 
deg + (p). Moreover, if p G then m(p) = m(p) = deg,(p). 

2.4. Proposition. Let p = a d x d + a<i-ix d ~ l + • • • + a G T> + (R) be of 
degree d > 2 and (r n ) nG No CE k a sequence given by a recurrence relation 
a-dfn+d + o d _ir n+d _! H h a r n = for n G No and r = ■ ■ ■ = r d _ 2 = 0, 

Td-l = 1- 

Then m{jp) = min{n G {d, d + 1, . . . }\r n < 0} 

We show now that nT(p) and n^(p) determine the value of deg + (p) for 
p G D+(R) of degree 2 or 3. 

2.5. Lemma. Let < ip < tt/2, < s < 1 and n = ^ . Piii = 
(cos kip, sin fcy>, (— s) fc ) G M 3 /or fc G No- T/iera 

(vfc x vo) • v/ = sin(/c — l)tp + (— s) fc sin lip — (s) sin kip. 

Let H(u,v) = {w G M 3 |(u x v) • w > 0} denote the halfspace determined 
by the (linearly independent) vectors u,v G M 3 . Similarly, put H°(u,v) = 
{w G M 3 |(u x v) • w > 0}. Lei e 3 = (0,0,1) G K 3 . 

(i) v 2 ,...,v n _i G fT°(v ,vi). 

(ii) Lfnisodd, then H(v , Vi)nH(v n _i, v )niT(vi, e 3 ) C (v , . . . , v n _i) R+ 
(hi) Lfniseven, then H(y , Vi)nif(v n _ 2 , v )nil(v 1 , e 3 ) C (v , . . . , v n _i} R . 

Proof. The first equation is easy to verify. Further, for I = 2, ... ,n — 1 we 
have (vi x vo) • v/ = — sin(/ — l)ip — ssinlip — (—s) 1 simp < —ssinlip — 
(sin(2 — l)ip — sin tp) < 0, hence Vo, . . . , v ra _i G H°(vq, Vi). 

Let u = (cos ?/>, sin ^, r) G i7(vo,vi) be such that tt/2 < ip < it and 
r G M. Choose 0^w£ (v ,u) R+ n (vi,e 3 ) R+ . Then F(u,v ) = i2"(w,v ) 
and ff(vi,e 3 ) = iT(vi,w) and it is not difficult to show that H(vq,v\) n 
jy(u,v ) nil(vi,e 3 ) = ff(v ,vi) n JT(w,v ) nlf(vi,w) = (v ,vi,w) R+ c 
(vo,vi,u) R+ . Now, choosing either u = v n _i or u = v n _2 we get the first 
inclusion in appropriate cases (i) and (ii). □ 

2.6. Theorem. Let p G D + (M) be of degree 2 or 3. Then deg + (p) = 
max{n^(p), n^(p)}. 

Moreover, letp = {x + c){x — oj){x — uj) G M\x\, where uj = r (cos p + i sin ip), 

c, r > and ip G (0, tt/2). Put n = — , < e = np — tt < p and 
s = min{^, £}. Then 

(i) n < deg + (p) < n + 1, 

(ii) deg + (p) = n if and only if sin(ip — e) < s ■ sine + (— l) n+1 s n • simp. 

Proof. Using 12.21 and 12.31 we only need to prove the case of p = (x + 
c)(x — uS){x — uj), where u = r(cos<p + i sin p), c, r > and cp G (0, tt/2). 
Now, (i) follows immediately from 12.21 Further, with help of substitutions 
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x = t ■ x', t > and x = 1/x' we have that deg + (p) = deg + (q) and 
max{n^(p), n^(p)} = max{n^(g), n^(q)} for q = (x + s)(x — ^)(x — ^), where 
s = min{f , ^} < 1. Put v fc = v fc (g) for G N . 

By I2.5f i) we have ^ = n < n^(q). If deg + (g) = n, then deg + (g) = 

xv (q) by 12.31 and det(vo, vi, v n ) < by I2,5f i). On the other hand, if 
> det(vo, vi, v n ) = (vo x vi) • v„, then we need to show that — v n G 
(vo, . . . , v n _i) K+ . Let H(-,-) and e3 have the same meaning as in 12.51 
Clearly, -v n G iJ(v ,vi) nH(vi,e 3 ). 

If n is odd then (v n _i x vo) • (— v n ) = sin tp + s n_1 sin e — s n sin((/3 — e) > 
by[231 hence -v n G F(v n _i, v ). 

If n is even then (v n _2 x vo)-(— v n ) = sin 2ip+s n ~ 2 sine+s n sin(299— e) > 0. 
Thus -v n G F(v n _ 2 ,v ). 

Putting it all together we conclude with — v n G (vo, . . . , v n _i) R+ by 
I2.5( ii).(iii) and deg + (q) = n. The rest is now easy and the condition in 
(ii) is only an equivalent version of the inequality (vo x vi) • v n < 0. □ 

2.7. Remark. The lower estimation n^(p) provides no additional infor- 
mation (i.e. n^(p) = deg(p)) in the case when both the linear and con- 
stant terms have coefficients of the same (non-zero) signs. In particular, 
the equality in 12.61 does not hold for degree 4. For instance, for p = 
{x + 2){x + l/2)(x — u)(x — uJ), uj = cos if + isinip, is n^(p) = n^(p) = 4, by 

12^1 while deg + (p) > | . 

Apart of this, the result in 12.61 (together with I2.4p provides a simple 
way (by choosing particular halfspaces) how to determine the value of deg + 
and we ask therefore for a possible generalization for polynomials of higher 
degrees. 

For an estimation of deg, for a general polynomial is useful to investigate 
this value for powers of polynomials. Eaton [5] provided deg + in the case 
when a quadratic polynomial vanishes on a root of unity (see 12. 8p . Following 
basically the same idea we can easily get the estimation in 12.91 

2.8. Theorem. [5j Let n, k G N and p = (x - e in / n )(x - e ~' m l n ) G R[x]. 
Then deg + (p k ) = kdeg + (p) = kn. 

2.9. Theorem. Let p G M[x] \R + [x] be a monic quadratic polynomial, i.e. 
p = (x — uj){x — ZJ), where uj = r(cos (p + i sin tp), r > and p G (0, ^). Let 
2 < k G N. Then 



(1) —1 <deg + (p*)<* 



r 



In particular, kdeg + (p) — k + 1 < deg , (p ) < kdeg + (p). 

The next proposition shows that both the upper and lower estimations in 
can be reached. 
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2.10. Proposition. For every n > 3, k > 2 there are monic quadratic 
polynomials pi,p2 such that 

(i) deg + (pi) = n and deg + (p^) = kdeg + (p\) and 

(ii) deg + (p2) = n and deg + (j?2) = ^deg + (p2) —k + 1. 



Proof, (i) follows from [2781 For (ii) we first take pi = (x — e ln ^ n ^)(x — 
e -m/(n-i)y By |2.2l there is / = Y17=o a i- x% divisible by P2 such that a\ > for 



every % = 0, . . . , n + 1. Now g 



„n-l 



+ l) k 1 / is divisible by pf, deg(g) 



kn — k + 1 and every i-th coefficient of g is positive for i 
Using continuity we get (p G R such that 

and {gjp^)p\ ^ R + [a;], where P2 = (x — 
deg+(p2) = kdeg + (p 2 ) - k + 1. 



kn 

V 



0, • • .,deg(g). 
■ k(n - 1) + 1 



e 1¥? )(x - e -1 *). Hence, by 



□ 



2.11. Example. Even if the lower and upper bound in l2.9t fTj) for a given 
polynomial differ, it may happen that the upper bound is achieved. Taking 
p = x 2 — x + 3, we have 5 the lower- and 6 the upper bound for p 2 in l2.9HfT1) . 

But deg + (p 2 ) = 6, since (x + c)p 2 = x 5 + (c- 2)x 4 + (7- 2c)x 3 + (7c- 6)x 2 + 
(9 - 6c)x + 9c £ R + [x] for any c > 0. 



The previous examples show that knowing only 



fc~ 

<p 



and A; 



is gener- 



ally not enough for determining deg + (p fc ) in 12.91 Nevertheless, based upon 
numerical computing, we conjecture that a better estimation than l2.9t [T|) 

+ 1 for p a quadratic polynomial. 



< 



might be true, namely deg + (|/ 

Now we turn to the general case. Using formal power series, Turan [19J 
has shown, that for a polynomial p G R[a;] of degree d with zeros outside the 
angle-domain {0 ^ z G C||arc(z)| < (po < n/2} is 



(2) 



deg + (p) < - 



7T 



¥0 



+ 1 



In the next theorem we make a slight improvement of this result and provide 
a simpler proof. 



2.12. Theorem. Let p G £>+(R) \R+[x] and d = deg(p). Let <p be the least 
positive argument of all non-zero roots of p. Then 



deg + (p) < 





d 




" 7T " 




< 




( 








.2- 




9?0 





+ d. 



Proof. Let p = p\ . . .p r q\ . . . q s be a decomposition into irreducible polyno- 
mials in M.[x], where r > 1, s > and pi G R[x] \ R + [x], qj G R + [x] for i = 
l,...,r,j = l,...,s. ByE2]wehavedeg + (p) < (d - 2r) + £- =1 deg+fe) < 
(d-2r) + r|"£l < LiJCfsrl -2) + d, since r < L#J - □ 



2.13. Remark. We show that the estimation in 12.121 can be reached for 
some polynomials, while the bound ([2]) is bigger in these cases. 



THE LOWEST-DEGREE POLYNOMIALS WITH NON-NEGATIVE COEFFICIENTS 7 



Let k,n G N. Put Pl = ((x - e m / n ){x 
d\ = deg(pi). Then, by 12.81 1S deg + (j>i) = 



e -i?r/n 

kn 





G 


R[x 


, <P0 = 


di 


( 


TV 




2 









ir/n and 



is less than k(n + 1) = -y 



+ 1 



Similarly, for p 2 = pi(x + c), d 2 = deg(j?2) and c > big enough we have 



D(pi) = D(p 2 ) and hence deg + (p2) = kn + 1 
the estimation in ([2]) is (fc + ^)(n + 1) = - 



d 2 


( 


7T 


- 2 ) 


2 




<^0 





2 ) +d 2 , while 



+ 1 



If a polynomial is given only by its coefficients, it may be difficult to 
determine arguments of its zeros. Therefore it is natural to look for an 
estimation that will use other features of a polynomial. In a recent paper 
Zaimi [23] gave a bound of deg + (p) (based on estimation of Dubickas [21 E]), 
for p being irreducible over Q[x], using the discriminant and the Mahler 
measure: 

2.14. Theorem. [23, 1.1] Let p = \\^ = i{x — cti) be an irreducible polynomial 
in Q[x] with ^ a% G C \ M + for every % = 1, . . . , d. Then 

deg + (p) < 



arcsm 



M d ~ 



|A| \ 

d d+3 ) 



where A = Yli<i<j<d( a i~ a j) 2 ^ s ^ e discriminant ofp and M = Y\ i=1 max{l, jcKj]} 
is the Mahler measure of p. 

Following the idea of [23], we improve this result to get approximately 4 
times better bound. The next lemma is a corollary of a theorem of Dubickas 
®. 

d 

2.15. Lemma. [3] Letp = Y\ ( x —&i) G WL[x] be a polynomial with a%, . . . , a d G 

i=i 

C such that \a\\ > \a 2 \ > • • • > \a-d\- Let u G C \ M + be a root of p. Then 



\0J 



to 



1 

> — 
- T 



IAI 



2|w| " T V d d + 3 
where A is the discriminant of p and T = ^a^l 2 ■ ■ ■ |<ay— 1|- 

If (po is the least positive argument of all the roots of p from 12 .151 and ipo G 
(0,7r/2), then T < ([\i=i max{l, |ai|}) d_1 = M^ 1 and we have smtpo > 

M l-i \J JhtI by 12.151 The next corollary now follows by fl2 . 12f) . 

2.16. Corollary. Let p = Yli=o a i x% e ^ + \ be a monic squarefree 

polynomial, d = deg(p), p(0) ^ 0, A be the discriminant and M be the 
Mahler measure ofp. Moreover, let L be an upper bound for the magnitude 
of any root ofp. Then 

d 



deg + (p) < 



7T 



arcsm 



d d+3 J 



2)+d 
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forKe {M, L d l 2 }. 

As a bound in 12.161 may serve for instance L = 1 + max{|o n _i|, . . . , |oq|} 

Of course, the previous assertion depends heavily on the stage of knowl- 
edge about estimation of the distance of roots of a polynomial. Moreover, 
it may happen that the bound in 12.151 will be depressed by a small angle 
between other roots that do not influence the value of deg + (p). However, so 
far 12.161 is the best known criterion. 

3. An unusual property 

In the end let us introduce a conjecture that was inspired by an approach 
using formal power series. For a formal power serie P = X^£o aiX% ^ ^[M] 
and n £ No denote (P) n = XT=o a i xl e ^[%]- As a consequence of a proof 
in [19] we get that for an irreducible p G R[x] of deg ree2isp-(|) n G K+[x] 
for n = deg + (p) — 2 = D(p). 

On the other hand, for q = x 2 - x + \ is deg + (q 2 ) = 8, D(q 2 ) = 4 by 121)1 
and q 2 -(ji)4 <£ K+[x] since jj = 4+16x + 32x 2 +32x 3 -16x 4 + - • • . However, 
we still have q 2 • (|) 2 • (|)2 G M + [x], since q ■ (|) 2 G R + [x\. With support 
of similar computational experiments we formulate a conjecture. We say 
that a polynomial p G such that p(0) ^ has a property (*) if there 

are p 1 , . . . ,p k G R[x] such that p = p 1 ■ ■ ■ p k , p { ■ (^-)d( Pi ) G M + [x] for every 

i = l,...,k and Z>(p) = D(pi). 

Conjecture: Let p G be a quadratic polynomial and p(0) ^ 0. 

Then p k has the property (*) for every k G N. 

3.1. Example. This example shows a cubic polynomial p = (x + ^)(x 2 — 
|x + 1) does not have the property (*). Assume the contrary. By 12.21 is 
deg + (g) = 5 for q = x 2 — \x + 1 and thus we have deg + (p) = 5, since 
(x 2 + x + \ )p = x 5 + i|x + § G R+[x]. Hence D(p) = 2 / + 3 = 
D(x + i) + -D(g). Therefore must be p ■ (|)d(p) GM + [x] by assumption. But 
i = 2-x + fx 2 - f x 3 + --- andp- (i) 2 = fx 5 - ^-x 4 + f x 3 + 1 £ M+[x], 
a contradiction. 
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